SELF-SIMILAR BLOW-UP IN PARABOLIC EQUATIONS OF 

MONGE-AMPERE TYPE 



C.J. BUDD AND V.A. GALAKTIONOV 

Abstract. We use techniques from reaction-diffusion theory to study the blow-up and ex- 
istence of solutions of the parabolic Monge-Ampere equation with power source, with the 
following basic 2D model 

(0.1) ut^ '\D^u\ + \u\P-^u in x R+, 

where in two-dimensions |Z?^u| = u^xUyy — (uxyY' and p > 1 is a fixed exponent. For a class 
of "dominated concave" and compactly supported radial initial data wo(x) > 0, the Cauchy 
problem is shown to be locally well-posed and to exhibit finite time blow-up that is described 
by similarity solutions. For p G (1,2], similarity solutions, containing domains of concavity 
and convexity, are shown to be compactly supported and correspond to surfaces with flat sides 
that persist until the blow-up time. The case p > 2 leads to single-point blow-up. Numerical 
computations of blow-up solutions without radial symmetry are also presented. 
The parabolic analogy of IjO.ip in 3D for which \D'^u\ is a cubic operator is 

ut = \D'^u\ -t- \u\P~'^u in R^ x R+, 

and is shown to admit a wider set of (oscillatory) self-similar blow-up patterns. Regional 
self-similar blow-up in a cubic radial model related to the fourth-order M-A equation of the 
type 

Ut = -\D'^u\ + in R2xR+, 
where the cubic operator \D^u\ is the catalecticant 3x3 determinant, is also briefly discussed. 

This is an earlier extended version of 6^, where, in particular, we present a survey on various 
M-A models; see Appendix A. 



1. Introduction: our basic parabolic Monge-Ampere equations with blow-up 

1.1. Outline. Fully nonlinear parabolic partial differential equations with spatial Monge- 
Ampere (M-A) operators arise in many problems related to optimal transport and geometric 
ffows [I7j, image registration [15], adaptive mesh generation [221IH], the evolution of vorticity 
in meteorological systems [T7|, the semi-geostrophic equations of meteorology, as well as being 
extensively studied in the analysis literature; see Taylor [SH Ch. 14,15], Gilbarg-Trudinger [101 
Ch. 17], Gutierrez [1^, and Trudinger-Wang [86], as a most recent reference. To describe such 
equations, we consider a given function u G C^(]R^), for which D^u denotes the corresponding 
N X N Hessian matrix D'^u = llwr t II, so that in two- dimensions, d = 2, 



;i.l) \D'^u\ = det D^u = u^xUyy - {u^y) 
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Similarly, in three dimensions 



(1.2) \D\ 

A general parabolic Monge-Ampere (M-A) equation with a nonlinear source term, then takes 
the form 

(1.3) ut = g{detD^u) +h{x,u,D^u) in x R+ 

with proper initial data u{0,x) = Uq{x). Such PDEs with various nonlinear operators g{-) and 
h{-) have a number of important applications. 

The origin of such fully nonlinear M-A equations dates back to Monge's paper [7^ in 1781, 
in which Monge proposed a civil-engineering problem of moving a mass of earth from one 
configuration to another in the most economical way. This problem has been further studied 
by Appel [Ij and L.V. Kantorovich [53l see references and a survey in pi]. Other key 
problems and M-A applications include: logarithmic Gauss and Hessian curvature flows, the 
Minkowski problem (1897) and the Weyl problem (with Calabi's related conjecture in complex 
geometry), etc. 

For increasing functions g{s), the equation fll.3p is parabolic if D'^u{-,t) remains positive 
definite for t > 0, assuming that V^uq > and local-in-time solutions exist [BT, p. 320]. 
Provided that g{s) does not grow too rapidly, for example if 

g{s) = Ins, g{s) = — and g{s) = for s > 0, 
it is known [STl HH [13] that the solutions of M-A exist for all time. 

In general, however, the questions of local solvability and regularity for M-A equations even in 
2D such as the system 

(1.4) {uxx + a){uyy + c) - (u^y + by = f 

in the hyperbolic (/ < 0) and mixed type (/ of changing sign) are difficult, and there are some 
counterexamples concerning these basic theoretical problems and concepts; see |55]. Note that 
classification problem for the M-A equations such as (11.41) on finding their simplest form was 
already posed and partially solved by Sophus Lie in 1872-74 ^1]; see details and recent results 
in [57]. 

For other functions g{s) in ( 11. 3p with a faster growth as s ^ oo and for certain nonlinear 
source terms h, solutions which are locally regular may evolve to blow-up in a finite time T. 
This gives special singular asymptotic patterns, which can also be of interest in some geometric 
applications; on singular patterns for M-A fiows, see [TSl [TU] The analysis literature currently 
suffers from a lack of understanding about the formation of such singularities in the fully non- 
linear M-A equation despite their relevance to such problems as front formation in meteorology 
[T7] . This paper aims to make a start at studying such blow-up singular behaviour by using 
techniques derived from studying react ion- diffusion equations to look at some special parabolic 
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M-A problems, which lead to the finite time formation of singularities. In particular, we con- 
sider f ll.Sp . fll.2p . and some other higher-order PDEs as formal basic equations demonstrating 
that M-A models can exhibit several common features of blow-up, which have been previously 
observed in PDEs with classic reaction-diffusion, porous medium, and the p-Laplacian oper- 
ators. Indeed, we will approach the study of M-A type operators by developing the related 
theory of the p-Laplacian operator. Our interest in this paper will be an understanding of the 
various forms of singularity that can arise in some parabolic Monge-Ampere models with a 
polynomial source term as well as extending the general existence theory for such problems. 

Further discussion of a variety of M-A parabolic, elliptic, and hyperbolic models together with 
basic regularity and singularity results continues in Appendix A after the list of references. 

1.2. Model equations and results. Model 1 Our first model fully nonlinear PDE is given 
by 

(1.5) Ut = {-lY-^\D\\ + \u\''-\ in R'^ x M+, 

where p > 1 is a given constant. Such models are natural counterparts of the porous medium 
equation with reaction/absorption, and of thin film (or Cahn-Hilliard-type, n = 0) models, 

(1.6) Ut = Au"'±uP and Ut = -W ■ {\u\''V Au) ± A\u\p-^u. 

Our principle interest lies in the study of those solutions which have large isolated spatial 
maxima tending towards singularities forming in the finite time T. Such solutions are locally 
concave close to the peak. The choice of sign of the principal operator in (11. 5p ensures local 
well-posedness (local parabolicity) of the partial differential equation in such neighborhoods. 
Significantly, the existence theory for such locally concave solutions is rather different from 
the usual theory of the M-A operator, which is restricted to globally convex solutions, and we 
will look at it detail in Section 4. The initial data Uq{x) > is assumed to be bell-shaped 
(this preserves "dominated concavity" ) and sometimes compactly supported. We firstly study 
radially symmetric solutions in two and three spatial dimensions, and will show analytically, by 
extending the theory of p-Laplacian operators, and demonstrate numerically, that whilst the 
Cauchy problem is locally well-posed, and admits a unique radially symmetric weak solution, 
certain of these solutions become singular with finite-time blow-up. We will also find a set of 
self-similar blow-up patterns, corresponding to single-point blow-up if p > d, regional blow-up 
for p = d, and to global blow-up for p G {l,d). The stability of these will be investigated 
numerically, and we will find that monotone self-similar blow-up profiles appear to be globally 
stable. 

We will also present some analytic and numerical results for the time evolution of non radially 
symmetric solutions in two dimensions. These computations will give some evidence to conclude 
that stable non-radially symmetric blow-up solution profiles also exist, though this leads to a 
number of difficult open mathematical problems. 

Model 2 As a second model equation, we will look at higher-order fully nonlinear M-A spatial 
operators associated with the equations of the form 

(1.7) ut = i-iy~'^\D\\ + \ufu in M'^ x M+, 
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where \D^u\ is the determinant of the 4th derivative matrix of u. We will find similar results on 
the blow-up profiles to those for the second-order operator. A principal feature of compactly 
supported solutions to (11. 7p is that these are infinitely oscillatory and changing sign at the 
interfaces, and this property persists until blow-up time. 

The layout of the remainder of this paper is as follows. In Section 2 we look at single-point blow- 
up, regional blow-up, and the global one of the radially symmetric solutions of the polynomial 
M-A equation in two and three spatial dimensions. We will combine both an analytical and a 
numerical study to determine the form and stability of the self-similar blow-up solutions. In 
Section 3, we will extend this analysis to look at solutions, which do not have radial symmetry, 
and will give numerical evidence for the existence of stable blow-up profiles in non radial 
geometries. In Section 4, we will look at more general properties associated with M-A type 
flows, in particular the existence of various conservation laws. Finally, in Section 5 we will 
study the forms of the blow-up behaviour for the equations with higher-order operators as in 

2. Parabolic M-A equations in M^: blow-up in radial geometry 

2.1. Radially symmetric solutions: first results on blow-up. The Hessian operator 
restricted to radially symmetric solutions in takes the form of a non-autonomous version 
of the p-Laplacian operator. Namely, for solutions u = u{r,t), with the single spatial variable 
r = \{x,y, ...)\ > 0, equation (11.51) takes the form 

(2.1) Ut = i-iy-'^ {Ur)''-'Urr + \u\P-'u in R+ X R+, 

and then for r = we have the symmetry condition 

Ur = 0. 

In this section, we shall consider the nature of the blow-up solutions for this problem and 
will identify different classes (single-point, regional and global) of self-similar radial solutions, 
giving some numerical evidence for their stability. However, we note at this stage (and will 
establish in the next section) that (possibly stable) non-radially symmetric blow-up solutions 
of the underlying PDE also exist. One can see that (12.11) implies that the equation is (at least, 
degenerate) parabolic if 

, , 1 , , , 1 f Mr < for even d, 

2.2 {-lY-\urY-'>0 ^ \ r 

[Ur IS arbitrary for odd a. 

For the local well-posedness of the above M-A flow, (12. 2p is always assumed. In all the cases, 
the differential operator {—lY~^ {urY~^Urr is regular in the class of monotone decreasing, 
sufficiently smooth, and strictly concave at the origin functions, so that the local well-posedness 
of (12.11) is guaranteed for the initial data satisfies the regularity and monotonicity constraints 

(2.3) u{r, 0) = Uo{r) > 0, uq E C\[0, oo)), Uq{0) = 0, M'o(r) < for r > 0. 
The corresponding p-Laplacian counterpart of (12. ip is then 

(2.4) ut = ^ lurl'^-^Urr + \u\P~^u in M+ X M+, 



which is locally well-posed by parabolic regularity theory; see e.g., p3]- By the Maximum 
Principle (MP), the assumptions in (12.31) guarantee that the solution u = u{r,t) satisfies the 
monotonicity condition 

(2.5) Ur{r,t)<0. 

Therefore, equations (12.11) and (12. 4p coincide in this class of monotone solutions. Note that for 
local well-posedness, we do not need any concavity-type assumptions that are usual for standard 
M-A fiows. 

The phenomenon of blow-up for the solutions of the model (12. 4p can be studied by using 
techniques derived from the study of reaction diffusion equations (see [HI Ch. 4]). By a 
comparison of the solution with sub-and super- solutions of the same equation of self-similar 
form, we can show that, for nonnegative solutions u, there exists a critical Fujita exponent 

(2.6) po = d + 2 such that: 

(i) for p G {l,Po], any u{x,t) ^ blows up in finite time, and 

(ii) in the supercritical range p > po = d + 2, solutions blow-up for large enough data, while 
for small ones, the solutions are global in time. 

The proof of blow-up in the critical case p = po is most delicate and demands a monotonic- 
ity /asjTuptotic rescaled construction; see e.g., [33l [36] . 

For the remainder of this paper we shall always assume that the initial data are such that 
blow-up always occurs. 

2.2. Blow-up similarity solutions. The M-A equation with the source term is in- 

variant under the scaling group 

p-d 1 

t ^ Xt, r A^'^Cp-i) r, u^X p-^u. 

Accordingly, a self-similar blow-up profile (with blow-up at the origin r = 0, which is assumed 
to belong to the blow-up set) is described by the following solutions: 

(2.7) Mr,t) = (T-t)-^/(z), z=(^, /3=2|^- 
Here / > is a solution of the following ordinary differential equation, 

(2.8) ^ (-i)rf-i(/')<^"i/" - ^fz - ^ / + l/rV = 0, f\z) < 0, /'(O) = /(+oo) = 0. 

The condition on /(+oo) (and the consequent requirement that the solutions of (12.81) should 
decay to zero as z ^ oo) is necessary to ensure that the self-similar solutions correspond 
to solutions of the original Cauchy problem of the PDE. In the case of monotone decreasing 
solutions with f'{z) < the equation (12. 8p becomes 

(2.9) \fV-'f" - Pf'^ - ^ / + Ifr'f = 0, /'(^) < in M+; /'(O) = /(+oo) = 0. 

In particular, in the case of (i = 2 to be studied in greater detail, we have to require the 
monotonicity assumption to allow the construction of smooth solutions. In the case of c? = 3, 
we can relax this assumption, leading to a richer class of (possibly oscillatory) solutions. When 
considered as an initial value problem, the solutions of (12.80 lose regularity at the degeneracy 
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/' = (except the origin r = 0) and are not twice differentiable at such points. However the 
existence of weak solutions with reduced regularity is guaranteed by the standard theory of the 
p-Laplacian operator, and these questions are standard in parabolic theory; see [23] and [391 
Ch. 2]. Note that the ODE (12.91) has two constant equilibria given by 

(2.10) ±/o = ±(p-l)-^^ 

and that solutions close to these equilibria can be oscillatory, which is a crucial property to be 
properly treated and used. 

If p > 1, then each solution of the form (12. 7p blows up in finite time, however the nature of 
the scaling is very different in the three cases of l<p< d, p = d and p > d, corresponding 
to global (blow-up over the whole of Mf^) , regional (blow-up over a sub-set of Mf^ with non-zero 
measure), and single-point blow-up respectively (zero measure blow-up set in general). Indeed, 
we will show that for p G {l,d] the solutions f{z) of the ODE are compactly supported, while 
for p > d they are strictly positive. 

2.3. Regional blow-up when p = d. We begin with the case p = d, where, according to 
(12. 7p . the ODE becomes autonomous and z = r, so that 

(2.11) ^ l/r^V" - / + l/rV = 0, fir) < in M+; /'(O) = /(+oo) = 0. 

This problem falls into the scope of the well-known blow-up analysis for quasilinear reaction- 
diffusion equations, [SH Ch. 4]. 

Proposition 2.1. The problem i \2.11\i has a non-trivial, monotone, compactly supported solu- 
tion Fo{z) > such that Fo(0) > 1. The support of Fq{z) is given by [0, Ls] with the asymptotic 
behaviour near the interface: as z —>■ Ls, 

(2.12) Fo{z) = AiLs-z)t'il + o{l)), where A = [^^^]—\ 

Proof. The result follows the lines of the ODE analysis in [STl pp. 183-189] and uses a shooting 
approach in which (12.81) is considered as an IVP with shooting parameter /(O). If /(O) is too 
large then the solutions of the IVP diverge to — oo and if it is not large enough then the solution 
has an oscillation about /o in a manner to be described in more detail below. The self-similar 
solution occurs at the point of transition between these two forms of behaviour. □ 

The form of the proof leads to a numerical method based on shooting for constructing an approx- 
imation to the function Fq. To do this we specify the value of Fq{0), take Fq{0) = and solve 
(12.111) as an initial value problem in r using an accurate numerical method (typically a variable 
step BDF method). In this numerical calculation the term |/'| is replaced by \/£^~+Xr)^ with 
e = 10~^. This allows the numerical method to cope with the loss of regularity when /' = 0. 
The value of /(O) is then steadily increased from /q and the transition between oscillatory and 
divergent behaviour determined by bisection. At this particular value /(O) = -Fo(O) the solution 
approaches zero as z Ls and for z > Ls v^e have Fs{0) = 0. Unlike studies of reaction- 
diffusion equations, the proof of the uniqueness of the solution of (12.111) is not straightforward. 
However, the numerical calculations strongly indicate the conjecture that such a compactly 
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Fj,{0)=1.814279.. 




supported monotone decreasing profile Fq{z) is indeed unique. In the case of d = 2, its support 
is 

(2.13) L5 = 3.26..., with Fo(0) = 1.814279... . 
Similarly, in the case of d = 3, 

(2.14) L5 = 2.303..., with Fo(0) = 1.366... . 

In Figured], we show the resulting profile (bold) in the case of p = d = 2, obtained by shooting 
as described above, together with an oscillatory solution of the IVP when /(O) ~ /o and some 
nearby divergent solutions. 

In the case of d = 3, we can relax the monotonicity requirement on the function f'{z). In 
this case there exists a countable set {F^} of compactly supported profiles that change sign 
precisely k times for any = 0, 1, 2, ... . A numerical shooting calculation of both Fo{z) (dotted) 
and of F[{z) (bold), together with some nearby oscillatory solutions, is shown in Figure [2l This 
figure both explains how these further profiles can be obtained numerically and indicates how 
their existence can be justified rigorously along the lines of the proofs in [5[ |35]. Here, we have 

if = 3.95... , with F[{0) = 1.6513... . 

It follows immediately that the variable separable solution given by 

(2.15) M5(r,t) = ^Fo(r) 

describes regional blow-up which is localized in the disc/ball {r < Ls}- 

We now make a further numerical calculation to investigate the stability of such a blow-up 
profile (in the restricted class of radially symmetric solutions). To do this we use a semi-discrete 
numerical method in which we discretise the Monge-Ampere spatial operator on a fine spatial 
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Two P-profiles for p=3, F^(0)=1.6513... 
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Figure 2. p = d = 3: the similarity profiles Fq (dotted) and Ff (bold) and some 
nearby oscillatory solutions, obtained by shooting in the ODE ()2.1ip from the origin 
with the parameter of shooting /(O) > 0. 

mesh. This leads to a set of ordinary differential equations for the discrete approximation to the 
solution u{r,t). These (stiff ordinary differential) equations are then solved using an accurate 
variable order BDF method. In this calculation we substitute a spatial domain r G [0,L] for 
the infinite interval and impose a Neumann condition at the boundaries r = and r = L. For 
a calculation with p = d = 2 we take L = 8 and use a spatial discretisation step size of L/1000. 
We present in the following figures some calculation showing the evolution of u{r,t) and the 
scaled function u{r,t)/u{0,t) taking as initial data 

u{r, 0) = 10e-°"'/l 

We consider two values of a to give profiles which lie above and below the self-similar solution. 
With a = 0.1 the solutions initially lie above the self-similar solution and in this case we see 
clear evidence in Figures [3] and H] for evolution towards self-similar regional blow-up with a 
blow-up time of T ^ 0.1099. 

We also plot in Figure O the value of m(t) = l/u(0,t). For small values of m this figure is very 
close to linear, and a linear fit gives 

m{t) ^ -0.5553t + 0.0610, so that u{0, t) ^ 

which is in good agreement with the earlier calculation of the self-similar profile. 
For comparison we now take initial data a = 10. In this case the value of u{0, t) initially 
decreases, and then increases with a blow-up time of T ^ 1.4371. Asymptotically the blow- 
up is almost identical to that observed earlier, however we can see clearly that this time the 
function u{r,t)/u{0,t) approaches the regional self-similar profile from below. 
Very similar figures arise in the case of d = 3 indicating that the regional self-similar solution 
is also stable in this case. 
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Figure 3. p = d = 2: Regional blow-up of the function u{r, t) with a = 0.1. 




Figure 4. p = d = 2: Regional blow-up of the scaled function u{r,t)/u{0,t) with 
a = 0.1 showing convergence to the self-similar profile with support [0,3.26] from 
above. 




Figure 5. p = d = 2: Evolution of l/n(0,t) with a = 0.1 showing that u{0,t) 
increases towards infinity. Note the linear behaviour ^ T — t close to the blow-up time. 



Figure 6. p = d = 2: Regional blow-up of the scaled function u{r,t)/u{0,t) with 
a = 10 showing convergence to the self-similar profile with support [0, 3.26] from below. 




Figure 7. p = d = 2: Evolution of 1/'U,(0, t) with a = 10 showing that ^(0, t) initially 
decreases before tending towards infinity. Note the linear behaviour ^ T — t close to 
the blow-up time. 

2.4. Single point blow-up for p > d: V and Q profiles. \i p > d then /5 > and single 
point blow-up occurs at the origin. The self-similar blow-up profiles can then take various 
forms. Initially we consider the monotone profiles for general d. 

Proposition 2.2. For any p > d, the ordinary differential equation problem i\2.S\) admits 
two strictly positive solutions Fff{z) > and Fq{z) > 0, which each satisfy the asymptotic 
expansion 





Fo(0) > /o, and 
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p=3, existence of a connection witti f=0; Fjj{0)=0.977513... 
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Figure 8. Similarity P-profile for d = 2 and p = 3 in (j2.9p obtained by shooting 
from the origin y = with the parameter of shooting Fq{0) > 0. 

(ii) Fq{z) = fo on some interval z G [0,ao] {i-e., it has flat sides in this bait), has the 
behaviour close to the interface given by 

(2.18) Foiz) = fo - ^ Paliz - ao)l{l + 0(1)) as z ao, 
and is strictly monotone decreasing and smooth for z > ao- 

The main ingredients of the proof of both P-type profiles of the form (i) and Q-type profiles 
of type (ii) are explained in [351 E]- Both solutions are again constructed by shooting, with 
/(O) being the shooting parameter for the P-type profiles and ao for the Q-type profiles. In 
Figure [H we show the similarity profile of the P-type solution Fo{z) (bold) for the case of 
d = 2 and p = 3 together with some nearby divergent solutions. In this case we find that 
Fo(0) = 0.9751... . Numerically there is strong evidence for uniqueness of this solution, but a 
proof of this uniqueness result is open. 

Similarly, in Figure [HI we show the results of shooting to find the Q-type profile when for = 2 
and p = 3 and we obtain numerically that ao = 2.292... . 

In the case of = 3, we may extend these results to construct countable families of non- 
monotone P and Q-type solutions described as follows: 

Proposition 2.3. If d = 3 then for any p > 3, problem {\2.8\] admits the following two countable 
families of solutions: 

(i) A P-type family {F[{z) > 0, A; > 0} such that 

2d 

(2.19) Ffc(O) ^ fo and Fk{z) = Coz"^(l + o(l)) as z ^ +oo, 

where Ck > is a constant, and each Fk{z) has precisely k + 1 intersections with the constant 
solution fo; and 
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Figure 9. Similarity Q-profile for d = 2, p = 3 in (|2.9p obtained by shooting from 
z = oq > 0, which is the parameter of shooting. 

(ii) A Q-type family {F^{z), k > 0}, where each F^{z) = fo on some interval z G [0,afc], 
{ttfe > 0} is strictly monotone decreasing, with the following behaviour at the interface: 

(2.20) F,(z) = /o + yp^(^-afc)|(l + o(l)) as z ^ a,,, 

and has precisely k intersections with /q for z > and has the asymptotic behaviour {\2.19^ . 
For the main concepts of the proof, see [351 E] • 

In Figure [TU] (a,b), we show the similarity profiles Fq{z) F^i^z) for d = 3,p = 4. Construction 
of further P-type profiles is similar and the following holds: 

F,^(0)^/o as A;-.oo, 

and moreover the convergence is from above for even k = 0,2,4, and from below for odd 
k = 1, 3, 5, ... . Uniqueness of each Fk with k + 1 intersections with equilibrium fo is a difficult 
open problem. We claim that as — > oo, both families {F^{z)} and {F^{z)} satisfying 

Ff(0)^/o and a? 
converge to a unique S-type profile F^{z) such that 

(2.21) F^(0) = fo and has infinitely many intersections with /o for small 2; > 0. 

The first two similarity profiles of Q-type with d = 3 and p = 4, with a fiat centre part, are 
shown in Figure [111 

Passing to the limit t T~ in (12. 7p . where 2 — > +00, for any fixed r > 0, we find from (12.171) 
the following final-time profile for both P and Q patterns: 

(2.22) us{r,T ) = Cor < 00 for any r > 0, 
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P-profile for p=4, F^(0)=1 .078739315. 



P-profile for p=4, F^(0)=0.4553... 




2 2.5 3 3.5 4 1 .5 2 2.5 3 3.5 4 

z z 

(a) /fc = (b) fc = 1 

Figure 1 1 . Q-type profiles of the ODE ^Mj for p = 4. 

with Co > fixed in (12.161) . This imphes single point blow-up at the origin r = only in both 
cases. 

We now make a numerical study of the stability of these solutions with single-point blow-up. 
Using a similar method to that described in the previous section (including taking Neumann 
boundary conditions with L = 8) we can study the nature of the time dependent blow-up 
solutions in this case. Usually when studying single-point blow-up the narrowing of the solution 
peak as the reaction time T is approached, requires the use of a spatially adaptive mesh |7j. 
However in the M-A systems, as the blow-us the width of the solution peak scales relatively 
slowly (as (T — tY^^ when d = 2) it is not necessary to use an adaptive method to study the 
nature of the blow up solutions in this case provided that the spatial grid is fine enough. Taking 
d = 2 and p = 3 and initial data u{r, t) = 10e~^ we find that the solution blows up in a time 
T = 0.028476. Indeed, plotting 1/^^(0, t) as a function of t we find a close to linear solution, 
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Figure 12. The solution u{r, t) when d = 2 and p = 3 as it evolves toward a P-type 
singular solution at time t = 0.028476 
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Figure 13. The rescaled solution (T — t)^/^u(r, i) plotted as a function of the sim- 
ilarity variable y = r/{T — t)^^^ compared with the similarity profile of the P-type 
solution Fq obtained by shooting. 

which has a best fit with the equation ^(O, t) = 0.975/(T — t)^/^ with T as above. Plotting the 
rescaled solution (T — ty^'^u{r,t) as a function of the similarity variable z = r/(T — t)^/® we 
find close agreement to the similarity solution constructed above. This strongly implies that 
the monotone P-type similarity solution is stable in the rescaled variables. A similar result is 
observed for calculations when d = 3. However, the Q-type and non-monotone P-type blow-up 
profiles all appear from these calculations to be unstable. 

2.5. S-type periodic solutions. The main part in the existence analysis for (12.91) . following 
the methods described in [351 [5], relies on constructing two solutions to the initial value problem 
with different oscillatory structure, and then deducing the existence of an intermediate solution 
with the correct properties at infinity by applying continuity arguments. To construct such 
solutions we must determine the oscillatory properties of the solutions f{z) about the positive 
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equilibrium /q. To study these we consider solutions of the form 

(2.23) f{z) = fo + yiz), where Y is small and solves the reduced equation 



(2.24) j^\Y'\'^-^Y" - pY'z + Y = 0, 

which remains nonlinear. In order to study oscillations of Y{z) about 0, we can exploit the 
scaling invariance of (12.241) and introduce the oscillatory component ip as follows: 

(2.25) Y{z) = z'^^is), s = hiz, where a = j^. 

Substituting fl2.25p into fl2.24p in (for example) the case old = 2 yields the following autonomous 
ODE: 

(2.26) + + 1^' + 12^) - + ^ = 0. 
Setting ip' = ipi^f) yields the first-order ODE system 

^^■''^ ^d^ = 1^ + 4^1 - - '2c,' 

which we can study by using phase-plane analysis, and particular identify a limit cycle. 



Proposition 2.4. The ODE l \2.27 ] admits a stable limit cycle on the {ip,ip} -plane, which 



generates a periodic solution 'p^:{s) of i\2.26^ . 

Proof. Equations (12.271) from PME and p-Laplacian theory with limit cycles have been studied 
since 1980's; see [35] and extra references and related results in [5j. The limit cycle exists for 
all p> d. □ 

As an immediate consequence, the gradient-dependent ODE (12. 9p also admits an S-type solution 
F^{z) with infinitely many oscillations about the equilibrium /o near the origin. Since the 
behaviour (12.230 . (I2.25P violates the monotonicity, this S blow-up profile is not associated with 
the original M-A equation when d = 2 but does correspond to a possible solution when d = 3. 
The existence of then follows by shooting from the origin using the ID asymptotic bundle 
constructed using (12.230 and the representation (12.250 with the periodic solution ip*{s), i.e., 

(2.28) f{z) = fo + z^^ip4so + \nz) + .... 

Here the shift Sq is the only shooting parameter. Using this parameter, it follows from the 
oscillatory structure of the expansion (13.140 that there exists an sq such that / = F''^ has the 
power decay at infinity given by (I2.17P with a positive constant Co- We do not have clear 
evidence for the uniqueness of such F^{z) > and it appears from the numerical calculations 
that the associated self-similar blow-up profiles are unstable. 

2.6. Global blow-up for p G {l,d). Many of the mathematical features of the analysis of 
the similarity blow-up structures in this case are the same as for p > d, though the evolution 
properties are quite different. It follows from (12. 7p that p G {l,d) corresponds to the global 
blow-up, where 

(2.29) us{x,t) — > oo as t ^ T~ uniformly on bounded intervals in M+. 
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Figure 14. Shooting compactly supported similarity profile Fq{z) (the bold line) for 
p = 1.5 in ([23]); Fo(0) = 11.5785... . 

As in [81 J pp. 183-189], we obtain existence of a similarity profile. Since the PDE (12. ip is 
non-autonomous in space, uniqueness remains an open problem, since the geometric Sturmian 
approach to uniqueness (see [10] for main results and references) does not apply. 

Proposition 2.5. For any p G (1,2), problem i \2.9^ admits a compactly supported solution 



The numerical shooting construction of Fq is shown in Figure by the bold line for d = 2 and 



It is immediate that the M-A equation (11. 5p admits non-radially symmetric blow-up solutions 
which do not become more symmetric as t — T^. These solutions can be obtained directly 
from a non-symmetric transformation of the radially symmetric blow-up solutions described in 
the previous section and we describe below. However, it is unclear from the analysis whether 
these solutions are stable or not. In this section we consider these solutions and make some 
numerical computations to infer their stability. 

3.1. Regional blow-up for p = d = 2: the existence of non-radially symmetric blow- 
up solutions. For p = d = 2, (II. 5p is 




3. Examples of blow-up in a non-radial geometry in 




This partial differential equation admits self-similar blow-up solutions (12. 7p . i.e.. 
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Here the function f[x, y) now satisfies the following "elliptic" M-A equation, with decay to zero 
at infinity: 

(3.3) A{f)^-\D'f\ + f-f = in R\ 

The radial compactly supported solution Fo(r) described in Proposition 12.11 also solves ( 13. 3p . 
The existence of non-radial smooth solutions of (13.31) now follows from an invariant group of 
scalings. Indeed, if /(x, y) is a solution of (13.31) then so is the function 

(3.4) faix,y) = f{^,ay) for any a > 0. 

Indeed any rotation of this function is also a solution. Therefore, by taking the radial profile 
Fo(?/) supported in [0, L5] as a solution of (12.111) . we can obtain the non-radially symmetric 
blow-up solution 



(3.5) us{x,y,t) = ^Fo(v/(f)2 + (ay)2), 
This blows up in the ellipsoidal localization domain given by 

(3.6) Ea = {{x,y): (f)' + (ay)^ < L| (a^l)}. 
Note that its area does not depend on a: 

measEa = vrL^ for any a > 0. 

Thus, M-A equations such as (13.11) do not support the (unconditional) symmetrization phenom- 
ena, found for many semilinear and quasilinear parabolic equations; see [391 P- 50] for references 
and basic results. In classic parabolic theory, results on symmetrization are well known and 
are connected with the moving plane method and Aleksandrov's Reflection Principle; see key 
references in ^0} Ch. 9] and [391 p. 51]. However, all these approaches are based on the Maxi- 
mum Principle that fails for M-A parabolic flows like (13.11) . We conjecture that (13.31) does not 
admit other non-symmetric solutions but have no proof of this result. 

3.2. On the linearized operator. Checking the stability properties of the non-radial solu- 
tions of (13. 3p . one can easily derive the linearized operator about the radial state -Fo('") 

Moreover, since A is potential in L^, this Frechet derivative is symmetric, so there is a hope 
to get a "proper" self-adjoint extension of the linear operator (13.71) . Unfortunately, we should 
recall that (13. 7p cannot be treated as elliptic in the domain of convexity of Fq. In addition, 
since -Fo(r) is compactly supported, the operator (13. 7p have singular coefficients at r = and 
will be inevitably defined in a complicated domain with possibly singular weights, which makes 
rather obscure using such operators in studying the angular stability or unstability of the radial 
profile -fo(r). In any case, it is convenient to note that, due to the symmetry (13.40 . the stability 
is neutral, i.e., 

(3.8) 3 A = 0, with the non-radial eigenfunction V'o = ^ ^"L=i ~ ~ ^ 7^ -^o('") 
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(we naturally assume that the eigenfunction belongs to the domain of the self-adjoint extension). 
In other words, the stability/unstability will depend on an appropriate and delicate centre 
manifold behaviour. More precisely, if we perform the linearization /(r) = Fq + ^(t) of the 
corresponding to (13. 3p non-stationary flow 

(3.9) /r = A(/) =^ Y^ = A'{Fo)Y-\D^Y\+Y^. 
Then the corresponding formal centre subspace behaviouiO by setting 

Y{t) = a{T)^o + w^, w^±V'o (||w^(r)|| < |a(r)| for r > l) 

leads, on projection (as in classic theory, we have to assume at the moment a certain com- 
pleteneee/closure of the orthonormal eigenfunction subset, which are very much questionable 
problems), to the equation 

(3.10) d = 7oa2 t:^1, where 70 = {-\D'^Fo\ + F^, ^0). 

Therefore, for 70 > 0, stability/instrability of such flows depend on how the sign of 70 is 
associated with the sign of the expansion coefficient ^(t). 

More careful checking by using equation (13.31) and the eigenfunction in (13.81) of changing sign 
shows that 

(3.11) 7o = (Fo,V'o) = 0, 

so that this centre subspace angular evolution according to (13.101) is formally absent at all. Of 
course, this is not surprising, since our functional setting assumes fixing the domain {r < Ls} 
of definition of the functions involved, and this clearly prevents any angular evolution on the 
centre subspace that demands changing this domain. 

In view of such a non-justifying formal linearized/invariant manifold analysis, we will next 
rely on also rather delicate numerical techniques to check angular stability of blow-up similarity 
profiles and solutions. 

3.3. Single point blow-up in non-radial geometry: similarity solutions for p > d = 2. 

We can use a similar method to study non-radially symmetric single-point blow-up profiles. 
The self-similar solution (12.91) . 

(3.12) us{x,y,t) = {T-t)-^^fi^,r^), ^ = x/ (T - t)^ , = y/ {T - , P = 
now leads to a more complicated elliptic M-A equation 

(3.13) _|DV|-/5V/-C-^/ + |/r7 = in {C = {i,r^Y). 

As before, the group of scalings (13. 4p leaves equation (13.131) invariant. Therefore, taking the 
radial solution Fq{z) from Proposition 12. 2[ we obtain the family 

(3.14) F,(e,r/)=Fo(^(J)2 + (ar^)2) 



Existence of a centre manifold by standard invariant manifold theory [68j is a very difRcult open problem, 
which seems hopeless. 

18 



of non-radial solutions of fl3.13p (together with all rotations of these). However, this set of solu- 
tions may not exhaust all the non-symmetric patterns. To see this, we consider the linearization 
(12.231) about the constant equilibrium /o which leads to a nonlinear M-A elliptic problem: 

(3.15) _|DV|-/5V/-y + /^-[/55/,,-(/5,)V/^(/?^ + ^^) + / = in 

This fully nonlinear PDE does not admit separation of variables. We can see from (I3.13P that 
if f{z) ^ as z oo sufficiently fast, the far- field behaviour is governed by the linear terms, 

(3.16) -/5(/ce + M)-^/ + - = 0. 
Solving this gives the following typical asymptotics (cf. (I2.17P ): 

(3.17) f{z) = C{^)\z\-^^+... {z={^,vr), 

where C{fi) > is an arbitrary smooth function on the unit circle {|/i| = 1} in M^. The constant 
function Cq{^) = Cq > gives the radially symmetric similarity profile as in Proposition 12.21 
Furthermore the vr-periodic function C(/i) given in the polar angle (p by 

generates the ellipsoidal solutions (13.14p . We conjecture that other solutions are possible with 
Ci{(p) having smaller periods |, ... . However, at present, the existence of such is unknown. 

3.4. Numerical computations of the non-radially symmetric time-dependent solu- 
tions. We now consider a numerical computation of the blow-up profiles when Q = [0, 1] x [0, 1] 
is the unit square (so that d = 2), and we took p = 3. It is convenient in this calculation to 
impose Dirichlet boundary conditions. The PDE is solved by using a semi-discrete method for 
which the square is divided into a uniform grid (typically a 100 x 100 mesh) and the spatial 
Monge-Ampere operator discretised in space by using a second-order 9 point stencil. The re- 
sulting time dependent ODE system is very nonlinear and an implicit solver is very inefficient. 
Accordingly it was solved using an explicit, adaptive Runge-Kutta method with a small toler- 
ance. The discretisation in space leads to certain chequer-board type instabilities and these 
are filtered out at each stage by using a suitable averaging spatial filter applied to the ODE 
system. For initial data satisfying the Dirichlet condition we took 

m(0, t) = 10'^e~'^''%in(7rx) sin(7r?/), where = a^(a; — 0.4)^ + ^ (y — 0.6)^, 

for which a = 2 and x, y were a set of coordinates rotated at an angle of |. This data were 
chosen to have an elliptical set of contours close to its peak. 

This system led to blow-up of the discrete system in a computed finite time T ^ 1.17 x 10^^. 
(Note that the computed blow-up time decreases when the spatial mesh is refined). In Figure fT5l 
(a,b) we present the initial solution and its contours and in Figure [12] (a,b) a solution much 
closer to the blow-up time (so that it is approximately 10 times larger than the initial profile). 
Note that the elliptical form of the contours has been preserved during the evolution giving 
some evidence for the stability of the elliptical blow-up patterns. We also give in Figure [T7| 
a plot of the solution a slightly later time. Note in this case evidence for an instability close 

^We recall that the M-A flow under consideration is supposed to have some natural instabilities in the areas, 
where the concavity of solutions is violated; these questions will be discussed. 

19 



0.2 0.4 0.6 0.8 1 

(a) Contours 




(b) Profile 

Figure 15. Initial solution profile and contours for d = 2, p = 4. 

to a point where the solution profile loses convexity. It is not clear at present whether this 
is a numerical or a true instability. Certainly all of the numerical methods used had extreme 
difficulty in computing a significant way into the blow-up evolution. 



4. On conservation of mass, comparison, existence-uniqueness, etc. 

Here we briefly discuss some related analytical theoretical aspects M-A equations such as 
(13. ip . which appeared before and have been already slightly discussed some times. We skip the 
quadratic reaction term m^, and concentrate on the fully nonlinear M-A operator, for which we 
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Figure 16. Profile and contours of the solution for d = 2,p = 4 closer to the blow-up time, 
consider the Cauchy problem 

(4.1) ut = -\D\\ = -u^xUyy + {u^yY in x I^+j 0) = uq{x) > in M^, 

with sufficiently smooth compactly supported initial data satisfying some extra necessary con- 
ditions (e.g., of having "dominated concavity"). Since compactly supported solutions are sup- 
posed to be involved, we cannot use the advanced techniques of classic nowadays theory of 
convex (concave) solutions; see [IHl [51], etc. Notice that, even for convex solutions, local 
regularity theory for the basic M-A equation 

detD^ = /(x) > in a convex bounded domain i7 C M 
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Figure 17. The solution at a slightly later time showing a possible instability at as 
point where the profile loses concavity. 

is rather involved with a number of open questions; see [H] as a guide. On the other hand, 
even for 2D stationary M-A equations of changing convexity-concavity (see (11. 4p ). there are 
counterexamples concerning local solvability and regularity, [55] . These questions of M-A theory 
have not been developed in detail stilH. The study of finite regularity solutions of the simplest 
degenerate (at 0) M-A equation in with radial homogeneous f{x), 

(4.2) detD^u=\x\'' in Bi 

has some surprises [20]; e.g. for a > 0, there exist a radial and a non-radial C^'"^ solutions, 
for a G (—2, 0) the only radial does (this is about a delicate study of a single point blow-up 
singularity for (14. 2p ): see also [80] for regularity of radial solutions for (14. 2 p with the right-hand 
side /(| |a;p,M, | |VmP). Equation (14. 2 p has the origin in Weyl's classic problem (1916). 

We do not plan to discuss such delicate questions somehow seriously here (especially for our 
problems having so strong degeneracy and even changing of type), and restrict to first auxiliary 
aspects of such singularity phenomena for (14. ip . 

4.1. Source-type similarity solutions. These are easy to construct for the radial equation 

(4.3) Ut = ~lUrUrr, SO that 

(4.4) u4r,t)=t-'sF{y), y = ^, where F (y) = j-^ {d^ - y% d>0. 



Yet, it is remarkable that the basic question of whether there exist any examples of local nonsolvability, 
has remained open for this well-studied class of equations", fSS', p. 665]. 
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4.2. Scaling group: non-symmetric solutions. Let us mention that, obviously, (14. ip ad- 
mits a variety of non-radial solutions. Indeed, the equation is invariant under the following 
group of scaling transformations: 

(4.5) u{x,y,t)^^u{^^,ll), a, 6, c^O. 

Therefore, (14. ip does not support the symmetrization phenomena that are typical for many 
nonlinear parabolic PDEs. 

Using also the time-translation, we obtain from (14. 4p the following 4-parametric family of 
exact solutions: 

(4.6) u^x, y,t) = ^iT + t)-l - cl (T + t)--. + y^)] l- 

4.3. No order-preserving semigroup in non-radial geometry. We recall that in the ra- 
dial geometry, the semigroup for the parabolic equation (14.30 is obviously order-preserving. It 
turns out that in the non-radial setting, this is not the case. 

Proposition 4.1. In general, sufficiently smooth solutions of i \4 ■ 1\) do not obey comparison. 

Proof. We take two exact solutions from (14. 6p : u{x, y, t) with a = b = c = T = l and the general 
solution u^{x,y,t), and show that the usual comparison is violated in this family of non-radial 
solutions. Comparing positions of the interfaces at the x- and y-axes and the maximum values 
at the origin yields for initial data at t = that 

(4.7) u^ix, y, 0) < uix, y, 1) if f^T^ < 1, < 1, ^gf T"! < 1. 
On the other hand, the comparison is violated for large t ^ 1 if 

(4.8) ^(r + t)l>ti, i.e., ^>1. 

It is easy to see that the system of four algebraic inequalities in (14. 7p and (14. 8 p has, e.g., the 
following solution: 

a = l, b=l = 2, d = 2, Te(^,i). □ 

4.4. Towards well-posedness. We present the following very formal speculations, which nev- 
ertheless show that the non-fully concave M-A flow (14. ip is better well-posed as it can be ex- 
pected. Actually, this was observed in a number of numerical experiments discussed above. 
Assume that, due to an essentially deformed spatial shape of the solution (say, by means of 
choosing special "ellipsoidal" initial data), we consider the unstable area that is characterized 
as follows: -C \uxxUyy\ and, e.g., Uyy > cq > 0, i.e., the flow 

(4.9) ut = -UyyUxx + ••• , u{x, y,t) > (cf. ut = -cqu^x + •••) 

is backward parabolic with respect to the spatial variable x. Then, let us assume that the 
positive solution u{x,y,t) is going to produce a blow-up singularity in finite time as t ^ T~, 
and, say, let it be a Dirac's delta of a positive measure (blow-up of derivatives is a different 
matter, which sometimes can be also treated). Of course, we do not mean precisely that in this 
fully nonlinear equation, but can expect that a certain such tendency as t moves towards T 
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can be observed, as the linear PDE in tlie braces in (14.91) suggests. Hence, if such a tendency 
of approaching a ~ 6{x — xq) in x is observed, then, obviously, at this unstable subset 

(4.10) Uxx < —Co < ^> Ut = {—Uxx)Uyy + ... bccomes well-posed parabolic in y 

(here we again assume that Uxy does not play a role at this stage). In other words, such a 
simple localized pointwise singularity ~ 5(x — xq) in both variables x and y is unlikely. This 
means that the PDE (14.11) can exhibit a certain "self-regularization" even in the case of not 
fully concave data. We are not aware of any rigorous mathematical justification of such a 
phenomenon, and will continue to discuss this subject below using other arguments. 

4.5. e-regularization: on formal extended semigroup. We propose to construct a unique 
proper solution of (14. ip as a limit of a family of smooth regularized solutions {we} satisfying 
the regularized fourth-order uniformly parabolic equation 

(4.11) Ue-. ut = -eA'^u - \D'^u\ in x M+, 

with the same initial data Uq. 

Global and even local solvability of this CP for (14. lip is a difficult open problem. Here, 
A(m) = \D'^u\ is a potential operator in L'^{M?) with the inner product denoted by (■, ■). The 
potential is given by 

1 

$(m) = / (m, A{pu)) dp = I / u\D'^u\. 



Hence, equation (14. lip admits two integral identities obtained by multiplication by u and Ut, 

'^iIu^ = -eJ{Au)^-Ju\DM, 

(4.12) T 

/ Jiutf = E{u{T)) - E{uo), E{u) = -I JiAuf - i / u\D'u\. 



In particular, writing the last identity as 

(4.13) JJ{u,r + I liAuf + I / u\D\\ = -E{uo) 



we see that for a uniform bound on > it is necessary to have the following "dominated 
concavity property" : for all t > 0, 

(4.14) l-D^Uel > in domains, where > is not small. 

To support this, as a formal illustration, we prove an "opposite nonexistence" result: 
Proposition 4.2. Assume that, for smooth enough uq, 

(4.15) / Uo\D'^Uo\ < 0. 
Then Us{-,t) for e -C 1 zs not bounded in for large t > 0. 
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Recall that, for the problem (14.111) . the conservation of the mass (i.e., an uniform estimate) 
is available; see Section 14. 7[ Note that (I4.15P is not true in the radial case for decreasing 
Mo(r) ^ 0, since by (14.31) . 

/ Uo\D\o\ = J rwo = -| / {u'q)^ > 0. 

Proof. Estimating from the second identity in (14.121) 

Ju\Dhi\ < -^J{AuY + 3[lJ{Auoy + lJuo\D^uo\] 
and substituting into the first one yields 

(4.16) '^iJu'> lJ{AuY-'fJ{Auo)'-Juo\D'uo\ > -'^ J Uo\D'uo\ > 
for sufficiently small e > 0. Hence, under the hypothesis (I4.15p . 

(4.17) Ju^{t)>l\Juo\D^uo\\t^+oo as t ^ oo. □ 

It seems that the divergence (I4.17P of {u^} actually means that the approximated solution 
u{x,t) is not global and must blow-up in finite time. This would imply global nonexistence of 
solution of (14. ip if (I4.15P violates the dominant concavity hypothesis (I4.14p at t = 0. However, 
we do not know whether reasonable data satisfying (14.150 actually exist. For instance, the 
standard profiles in (14.60 do not obey (I4.15P (since they correspond to uniformly bounded 
L^-solutions of (lO) ). 

In general, identities (14.120 cannot provide us with estimates that are sufficient for passing 
to the limit as e — 0, so extra difficult analysis is necessary. The main difficulty is that the 
Hessian potential $(m) is not definite in the present functional setting and the operator A{u) 
is not coercive in the class of not fully concave functions. To avoid such a difficulty, another 
uniformly parabolic e-regularization may be considered useful such as, e.g., 

(4.18) Ue-. ut = -eA[{l + u^)Au] - \D'^u\. 

Unfortunately, the first operator is not potential in L^, so deriving integral estimates become 
more tricky. On the other hand, using degenerate higher-order p-Laplacian operators such as 

(4.19) Us-. ut = -eA(|Au|PAn) - \D^u\ 

can be more efficient for p > 1. Here both operators are potential in L^(M^), and moreover, the 
p-Laplacian one is monotone, which can simplify derivation of necessary estimates; see Lions' 
classic book [651 Ch. 2]. Nevertheless, using various e-regularizations as in (14.110 . (14.180 . or 
(14.190 does not neglect the necessity of the difficult study of boundary layers as e — » 0. 

Thus, according to extended semigroup theory (see [Mi Ch. 7]), the proper solution of the 
CP (14. ip is given by the limit 

(4.20) u{x, t) = \im,_,o Ue{x,t). 

In general, as we have mentioned, existence (and hence uniqueness) of such limits assumes 
delicate studied of e-boundary layers which can occur in the singular limit e ^ 0. In particular, 
the uniqueness of such a proper solution would be guaranteed by the fact that the regularized 
sequence {us} does not exhibit 0(1) oscillations as e ^ 0, so (I4.2O0 does not have different 
particular limits along different subsequences {et} 0. Another important aspect is to show 
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that the proper solution does not depend on the character of the e-regularizations apphed. 
Such a strong uniqueness resuh is known for the second-order parabohc problems [34) Ch. 6,7] 
and is based on the MP. For higher-order PDEs, all such uniqueness problems are entirely open 
excluding, possibly, some very special kind of equations. We hope that the potential properties 
of the Hessian operator \D'^u\ and hence identities like (14.121) can help for passing to the limit 
in (14. lip or other regularized PDEs and to avoid studying in full generality difficult singular 
boundary layers. These questions remain open. 

4.6. Riemann's problems: a unique solution via a formal asymptotic series. Let us 

continue to study the passage to the limit e — > in the regularized problem (14. lip . We assume 
that the origin belongs to the boundary of supp uq and 

(4.21) u^{X)=0{\\Xt) as X^O {X = {x,yf). 

This class of data specifies a kind of a Riemann's problem under consideration (with given type 
of singular transition to 0). Then we perform the following scaling in (14.110 : 

(4.22) u{X,t)=ev,{Ct), C = X/ei, 

so that fe(C,t) solves an e-independent equation with e-dependent data, 

(4.23) V,: vt = Ao{v) = -A'w - \D'w\, VosiO = 1 MC^'^) ■ 
According to (I4.2ip . we assume that there exists a finite limit on any compact subset 

(4.24) y^^^Q^y^^Q as e^O, 

where, without loss of generality, by fo(C) we mean a fourth-degree polynomial. 

As usual in asymptotic expansion theory (see e.g. Il'in [19]), the crucial is the first nonlinear 
step, where we find the first approximation Vo(C,^) satisfying the uniformly parabolic PDE 

(4.25) Vo: Vt = Ao{V), \/(C, 0) = t;o(C). 

It can be shown by classic parabolic theory [251 132] that the fourth-degree growth of fo(C) 
( ^ CO guarantees at least local existence and uniqueness of Vq. 

As the next step, we define the second term of approximation, v = Vq + wq, where wq solves 
the linearized problem 

(4.26) wq: wt = A'q{Vo)w, w{0,y) = voeiy) - voiy), etc., 

again checking that this linear parabolic problem is well-posed by classic theory [251 
Eventually, this means that we formally express the solution via the asymptotic series 

(4.27) v,{y,t) = Vo{y,t) + J2^^^Wj{y,t;e), 

where each term Wk for /c > 1, is obtained by linearization on the previous member, by setting 

Vk = Vk-i + Wk = Vo + J2j<k 
which gives for Wk a non-homogeneous linear parabolic problem 

(4.28) Wk-. wt = A'oiVk.i)w - [iVk-i)t - Ao(\4.„i)] , wiO, y) = 0, 
with similar assumptions on the well-posedness. 
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Thus, (14.271) gives a unique formal representation of the solution. In asymptotic expansion 
theory, the convergence of such series and passing to the limit e ^ are often extremely 
difficult even for lower-order PDEs, where the rate of convergence or asymptotics are also hardly 
understandable. For instance, the asymptotic expansion for the classic Burgers' equation 

(4.29) ut + uux = euxx 

contains Ine terms (close to shock waves), and a technically hard proof uses multiple reductions 
of fl4.29p to the linear heat equation, which is illusive for our M-A PDEs; see [49j. For practical 
reasons, it is important that, as an intrinsic feature of asymptotic series, each term in (14.271) 
(including the first and the simplest one wq) refiects the actual rate of convergence of Ue given 
by (14.201) as e to the proper solution (I4.20p . 

4.7. Conservation of mass. In view of (14.201) . this is straightforward, since 

(4.30) \D\\ = div V, where V = i ^"''^^ ~ ''^''"2/ 

Therefore, given initial data uq G L^(M), integrating (14. lip over M^, we obtain that 

(4.31) / U!r{x,t)dx= f Uo{x) dx for t > 0, 

R2 1R2 

meaning the conservation of mass of this regularized M-A fiow. In this sense, (14.111) looks like 
a fully nonlinear version of the quadratic PME from filtration theory 

Ut = Am^ = divgradw^, 

but corresponds to more nonlinear and complicated mathematics without standard MP and 
accompanying barrier techniques. 

5. Examples of blow-up for a fourth-order M-A equation with -\DSi\ 

5.1. On derivation of the higher-order radial M-A model. Our intention is to show 
that such models make sense and exhibit correct local parabolic well-posedness. Namely, we 
introduce radial models related to the fourth-order M-A equation 

(5.1) Ut = -\D%\ + u^ in x M, 

with the catalecticant determinant \D'^u\ given by 



(5.2) deW% = det 



'^xxxx 


'^xxxy 


U 


'^xxxy 


'U'xxyy 


U 


Uxxyy 


Uxyyy 


U 



^xxyy 
I'xyyy 

''yyyy _ 

which plays an important role in the theory of quartic forms. For instance, each such form in 
two variables can be expressed via a sum of three fourth powers of linear forms and via two 
powers, provided that det D'^u = 0; see 
It then follows from (15. 2p that 

1-^ ^1 ^xxxx^wvv^xxvv ~t~ '^^xxxv^xvvv^a 



(5.3) 



''xxxx^yyyy^xxyy ^ ^^xxxy^xyyy'-^xxyy 

(^xxj/j/) i,^xyyy) "^xxxx i^'^xxxy) ^yyyy 
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In particular, for radial functions u = u{r), we have 



6 I j,7 , 

y" ,,(4) I 6x2y2 ,„ 3(x4-4x2iy2) 3(-x4+4x2^^) , 



Uyyyy = ^ M^^^ + M'" + ^^^^^^^ M" + ;7'^ ^ W 

2, ,2 ^4_i_,,4_^^2,,2 „, _9^'^■l_l_„4^_l_11^2,,2 „ 9/'^4j_„4-,_i i ^2,,2 



(5 4) ^ a.'y'^ (4) I x*+^-^-4x^^2 „/ , -2(x^+i;4)+ll3;2^2 2(x4+y4) 11^2^- 

(4) I 3{a;y3-a;^y) /// , 3(2x^y-3xy^) // , 3(-2x'^y+3xy^) / 

7, _ xy^ U) I 3(x3y-xy3) ,„ 3(2xy3-3x3y) „ 3(-2xi;3+3x3y) , 

"■xyyy — j.4 " ~r ^5 u ~r ^6 " ~r ,,7 u • 

Overall, this leads to the following radial model that is the simplest in such flows: 

5.2. Simplest well-posed radial fourth-order M-A type model with regional blow- 
up. Using (15. 4p . by balancing and mutual cancellation of the most singular terms as r — 0, 
we introduce the following model radial parabolic equation associated with the M-A flow (15. ip : 

(5.5) Ut = —jj {UrrrYUrrrr + in IR+ X 

Looking at it as a parabolic equation, we pose at the origin the symmetry (regularity) conditions 

(5.6) UriS^^t) = Urrr{0,t) = 0. 

Obviously, (15. 5p is uniformly parabolic with smooth (analytic) solutions in any domain of non- 
degeneracy {urrr 7^ 0}. In particular, checking the regularity of the operator in (15.50 and 
passing to the limit r — yields 

(5.7) ^ ("^rrr) Urrrr ^ {,^rrrr) ; 

SO that, at the origin, the differential operator is non-degenerate and regular if Urrrr 7^ 0. Thus, 
regardless the degeneracy of the equation (15. 5p . this radial version of fourth-order M-A flows is 
well-posed, at least, locally in time. 

Describing blow-up patterns and looking for blow-up similarity solutions of (15. 5p . 

(5.8) Ws(r,t) = -^/(r), 
we obtain the ODE 

(5.9) A(/) = -^(r)¥4) + /3-i/ = for r>0, 
with the symmetry conditions generated by (15. 6p at the origin, 

(5.10) /'(O) = r (0) = 0. 

Indeed, (15. 9 P is a difficult fourth-order ODE with non-monotone, non-autonomous, and non- 
potential operators. Any (at least) 3D phase-plane analysis or shooting arguments are rather 
difficult. Instead, and this is much easier nowadays, we are going to check existence (and 
uniqueness) of the solutions using reliable numerical methods of an enhanced accuracy. For 
instance. Figure fTSl demonstrates the corresponding compactly supported solution of (15. 9p . We 
did not detect other P-type profiles that have more than one oscillation about the constant 
equilibrium 

•^0 = 72- 
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Figure 18. A unique compactly supported nonnegative solution of the ODE (15. 9p . 

As usual, stability of the similarity blow-up (15.81) is studied in terms of the rescaled solution 

(5.11) w{x,t) = VT-t u{r,t), r = - ln(T - t) ^ +00, 
where w solves the rescaled parabolic equation with the same elliptic operator 

(5.12) Wr = A{w) = {WrrryWrrrr + — in IR+ X M+. 

Since A is not potential in any suitable metric, so f l5.12p is not a gradient system, passage to the 
limit as r +oo to show stabilization to the stationary profile /(r) represents a difficult open 
problem. Note that the main operator — ^ {wrrrYwrrrr is gradient and admits multiplication 

by Wrrr in L'^i^), 

but (and also w) does not. 

5.3. p = 3: on oscillatory structure close to interfaces. In addition. Figure [18] shows 
that, locally, close to the finite interface point r = tq, sufficiently smooth solutions of (15.91) are 
oscillatory. This kind of non-standard behaviour of solutions of the Cauchy problem deserves 
a more detailed analysis. To describe this, we introduce an extra scaling by setting 

(5.13) f{r) = {ro — r)'^{p{s), s = ln(ro — r), where 7 = 3. 

Substituting (15.131) into (15. 9p and neglecting the higher-degree term for r ^ Tq , we obtain 
the following equation for the oscillatory component <f{s): 

(5.14) (f'yfi^) = -X,f =^ {PMfPM = -Xo^, Ao = ir2, 
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Figure 19. Stabilization to a stable periodic orbit (p^,{s) for the ODE (j5.14p . Aq = 1. 

where Pk are linear differential polynomials obtained by the recursion procedure (see p. 140]) 
Pi(y^) = / + -fifi, P2{ip) = if" + (27 - 1)/ + 7(7 - l)ip, 
Ps{y^) = + 3(7 - W + (372 - 67 + 2)y,' + 7(7 - 1)(7 - 2)^, 
Pi{^) = <^(^) + 2(27 - 3)^"' + (672 - I87 + 11)^" 
+ 2(273 - 972 + II7 - 3)ip' + 7(7 - 1)(7 - 2)(7 - 3)(^. 

In Figure [191 we show the typical behaviour of solutions of the second ODE in (15.141) demon- 
strating a fast stabilization to the unique stable periodic solution. According to (15.131) . this 
periodic orbit (f^{s) describes the generic character of oscillations of solutions of the ODE (15. 9p . 
Periodic solutions for semilinear ODEs of the type (I5.14p are already known for the third- [261 
§ 7] and fifth-order operators (with ^5(9?)) [271 § 6, 12]. The quasilinear equation (15.141) is more 
difficult, and existence and uniqueness of V5*(s) remain open. 

As a key application of the expansion (15.131) . we note that this shows that the ODE (15.91) 
generates a 2D asymptotic bundle of solutions close to interfaces 



(5.15) /(r) = (ro — r)'^y9(s + So) + ... as r 



'0 ) 



with two parameters, tq > and sq G M as the phase shift in the periodic orbit. The 2D 
bundle perfectly suits shooting also two boundary conditions (I5.10p . though a proper topology 
of shooting needs and deserves extra analysis. 

5.4. Structures of single point blow-up for p > 3. As usual, single point blow-up occurs 
for p > 3 in the M-A PDF with other power nonlinearities, 

(5.16) Ut = — {UrrrYUrrrT + iu M+ X M_j., 
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Figure 20. Single point blow-up profile satisfying the ODE ([5J81) for p G (3,8]. 

where we again pose the symmetry conditions (15. 6p . The blow-up similarity solutions are now 

(5.17) Ms(r,t) = (T-t)-^/(z), z = r/{T-tf, where /? = 
and / solves the ODE 

(5.18) A(/) = -^(r)¥4)-/?/'z-^/+|/r7 = for ^>0, 

with the same symmetry conditions at z = 0, ( ]5.10p . For p = 3, ( 15.18^ yields the simpler 
autonomous equation (15.91) for regional blow-up. For p > 3, this ODE is more difficult, and, 
following the results in Section [2], we expect to have blow-up profiles of P, Q and, possibly, 
S-type. 

Figure [201 demonstrates the first profiles of P-type for various p G [3,8]. The dotted line 
corresponds to the regional blow-up profiles, p = 3, from Figure [181 In particular, this clearly 
shows the continuity of the ( "homotopic" ) deformation of solutions of (I5.18P as p — 3^. 

It is key to observe that the profiles for p > 3 have infinite interface and, moreover, 

(5.19) f{z) > in M_|_ for p larger than, about, 5. 

For smaller p > 3, the profiles continue to change sign as for p = 3, as the continuity in p 
suggests. In Figure HH we show the enlarged behaviour of the profiles from Figure [201 in the 
domains, where these are sufficiently small. In both Figures (a) and (b), the profile for p = 3.5 
has two zeros only. 

It follows from (15.181) that the non-oscillatory profiles f{z) have the asymptotic behaviour 
governed by two linear terms (other nonlinear ones are negligible on such asymptotics), 

(5.20) _ _i_ / = =^ f{z) = Cz^ + ..., ^ = __i^<0, 
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Figure 21. Non-oscillatory behaviour of profiles from Figure [20] for p > 3. 

where C = C{p) for a.e. p > 3. The full 2D bundle of such non-oscillatory asymptotics 
includes an essentially "non-analytic" term of a typical centre manifold nature [73] : 

(5.21) /(z) = Cz^ + ... + Cie-'"'^^ + ..., where = > 0, 

Ci G M is the second parameter, and bo = bo{p) > is a constant that can be easily computed. 
The expansion (15.211) can be justified by rather technical application of standard fixed point 
theorems in a weighted spaces of continuous functions defined for large z ^ 1. 

The full 2D bundle ( ]5.2ip poses a well-balanced matching problem to satisfy two symmetry 
conditions at the origin (15.101) for any p > 3, provided that C ^ 0. 

On the other hand, Figure [2T](b) clearly shows that C (p) changes sign at some 

(5.22) PiG (4.5,5), and C{pi) = 0. 

Then, at p = pi, using the oscillatory analysis presented below, we expect that the correspond- 
ing f{z) is compactly supported. 

Moreover, in view of the oscillatory behaviour near interfaces (cf. (15.131) for p = 3), we expect 
that there exists a monotone decreasing sequence of such critical exponents 

(5.23) {pk}k>i ^ 3~^ as k ^ oo, such that C{pk) = 0, 

so that pi in (15.221) is just the first, maximal one. Possibly, such a mixture of compactly 
supported and non-compactly supported profiles f{z) for p > 3 gets more complicated (see 
further comments below on oscillatory character of finite- interface solutions). 

In the cases (I5.23p . the single point blow-up profile can be compactly supported, so we need 
to describe its local behaviour near the interface, which turns out to be different from that for 
p = 3 studied above. 

5.5. Finite interfaces: on oscillatory structures for p > 3. The oscillatory structure of 
solutions is given by two first terms in (I5.18p . and hence have a different form, than in (15.131) . 

(5.24) f{z) = {zq — zyip{s), s = ln{zQ — z), where 7 = |- 
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Figure 22. Stabilization to a stable periodic orbit ip^{s) for the ODE (|5.25p . 

Substituting into flS.lSp and neglecting other higher- degree terms, yields for ip{s) the ODE 

(5.25) (PMfP.icp) = -P,{cp) ^ -{cp' + 7y,), 

where we have scaled out the constant multiplier f3zQ > on the right-hand side. This ODE is 
even more difficult than fl5.14p . so we again rely on careful numerics. 

In Figure [221 "we show the typical behaviour of solutions of (15.251) demonstrating a fast 
stabilization to a unique and stable periodic solution. According to (15.241) . this periodic orbit 
V?*(s) describes the generic character of oscillations of solutions of the ODE (I5.18P in the critical 
case (I5:22|) . 

Since here the periodic orbit is stable as s — >• +oo, it is unstable as s — >• — oo (in the direction 
towards the interface at s = —oo according to (15.240 ). and moreover the stable manifold of 
(p^:{s) as s — oo consists of the solution itself up to shifting. Therefore, the full equation 
(I5.18P admits precisely 

(5.26) 2D bundle of small solutions, with parameters > and Sq G M, 

where Sq is again the translation in the oscillatory periodic component (f{s + Sq). Therefore, 
shooting via 2D bundle precisely two symmetry conditions at the origin (I5.10p represents a well- 
posed problem, which can admit solutions, and possibly a countable set of these. Nevertheless, 
Figures [2T](a) and (b) justify that the actual behaviour is governed by the non-compactly 
supported 2D bundle (I5.2ip for a.e. p > 3, and we do not know any mathematical reason why 
the finite-interface bundle (also 2D) fails to be applied for p > 3 a.e. and not only for p = pk. 

5.6. On multiplicity of solutions by branching. Figure [23] shows two P-type proffies, Fi{z) 
and F2{z). This poses a difficult open problem on multiplicity of solutions for p > 3 (and also 
for p = 3). Recall that operators in both equations (15.180 and (15. 9p are not potential, so we 
cannot rely on well developed theory of multiplicity for variational problems; see e.g., [59l § 57]. 
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Figure 23. Two P-type profiles of the ODE (f5l8]) for p = 4. 

Nevertheless, we will rely on variational theory by introducing a family of approximating 
operators 

(5.27) A^(/)= [-i^(n^ + /i-l]/W-/i/5A-^/+|/rV, /^e [0,1]. 
Then Ai = A, while for /i = 0, 

(5.28) Ao(/) = -/(^)-^/+|/r7, 

we obtain a standard non-coercive variational operator. The corresponding functional for / G 
H'^{—L, L) on a fixed interval with L ^ 1, 

(5-29) %{f) = -\ lif"? -^iP + ^^i \fr' 

has at least a countable set of different critical points {fj,^\ k = 0,1, 2, ...}; see 

Thus we arrive at the branching problem from profiles fj:^^ at the branching point /i = 0, 
which leads to classic branching theory; see [591 Ch. 6] and [88J . In general in the present ODE 
setting, where the linearized operator for /i = 0, 

(5.30) A'o(/)y = -r(^)+p|/riy, 

has a ID kernel, branching theory [2TI p. 401] suggests that each member fj!^^ generates a 

continuous branch {fjf"^} at /i = 0. The global continuation of those branches up to /i = 1 
remains a difficult open problem, that in the present ODE case admits an effective numerical 
treatment. 
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Appendix A: On various Monge-Ampere models and application 

Second-order Monge Ampere equations: critical exponents and singularities. Thus, for a 
given function u G C^(M^), D^u denotes the corresponding N x N Hessian matrix D'^u = \\uxix \\- 
As we have mentioned in Introduction, general parabolic Monge-Ampere (M-A) equations (jl.3p play 
a key role in many geometric problems and applications. Thus, M-A equations first appeared in 
Monge's paper [72J in 1781, related to civil-engineering problem of moving a mass of earth from one 
configuration to another in the most economical way. This was further studied by Appel \T\ in 1887 
and Kantorovich in the 1940s [53| I54j: see history in |31j . Other key problems and M-A applications 
are: 

(i) logarithmic Gauss and Hessian curvature flows, 

(ii) the Minkowski problem (1897), 

(iii) the Weyl problem (with Calabi's related conjecture in complex geometry), etc. 

We refer to basic mathematical results in monographs by Taylor [84^ Ch. 14,15], Gilbarg-Trudinger 
[40| Ch. 17], and Gutierrez ^J. For increasing functions g{s), equation (II. 3p is parabolic if D'^u{-,t) 
remains positively definite for t > 0, assuming that D^uo > for initial data uq. For a class of 
lower-order operators h{-) satisfying necessary growth estimates, typical monotone increasing, concave 
nonlinearities g in the principal operator are 

g{s) = Ins, g{s) = and g{s) = sn for s > 0, 

for which the global-in-time solvability is known. For other g{s) functions with a faster growth as 
s — > oo, local-in-time solutions existing by standard parabolic theory (see e.g., the classic book [61} 
p. 320]) may blow-up in finite time. This gives special asymptotic patterns, which can also be of interest 
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in some geometric applications. For instance, the affine normal flow for an initial convex, properly 
embedded, and noncompact hypersurface L C M^^^ can be describes by the parabolic equation 

(A.l) ut = -{detD^y-^; 

see [67] as a guide for the history and recent results including finite-time extinction. 
The parabolic M-A equation 

-Ut detD^u = f in Qt = 5^ X (0, T), 

where Q, C is a bounded smooth domain, was first introduced in [58j; see recent related references 
in [33] and more general models like that in [82]. Note the pioneering paper by Hamilton [l6] on the 
evolution of a metric in direction of its Ricci curvature. Conditions of the global unique solvability of 
the M-A equation 

Ut = (detD^ii)^ + g in Qt, 

were obtained by Ivochkina and Ladyzhenskaya [51]. This model corresponds to special curvature 
flows. 

The logarithmic Gauss curvature fiow in terms of the support function (see [13]) is described by the 
M-A-type equation 

Ut = Y^l + |x|2 ln{detD'^u) +h in x IR+, 

where h is given. Depending on the initial uniformly convex hyper surf aces, this PDE is known to 
admit either a local solution, corresponding to shrinking to a point in finite time or a global solution 
that describes the uniform convergence to an expanding sphere. 

The relations between Kahler-Ricci flows on M ((M",(7o) is a complete Kahler manifold) 

^ = -Rfj, gij{x,0) = {go)fj 

and the corresponding parabolic M-A equations 

dt - dot(9o)fcj) ■'0' ^ia;,Uj-U, 

is explained in [9], where further references on global solvability, asymptotics, convergence, and appli- 
cation of these problems can be found. 

Typical second-order Hessian operators are known to be potential and the corresponding smooth 
parabolic flows are gradient systems. For M-A PDEs, these ideas go back to Bernstein (1910) [3]; see 
also Reilly, and [74] for a particular case. Potential properties of more general Hessian operators are 
described in [87]. Parabolic M-A equations as gradient flows, and the related questions of the asymp- 
totic behavior of solutions, were studied, e.g., in [14:\ I82j. where further references concerning various 
types of Gaussian flows can be found. The classical Gauss curvature flow describes the deformation of 
a convex compact surface S : z = u{x,y,t) in by its Gauss curvature and is governed by the PDE 

(A.2) Ut = (YqUp^, 

which is uniformly parabolic on strictly convex solutions. Singularity formation phenomena for (IA.2P 
appear if the initial surface S has flat sides, where the curvature becomes zero and the equation degen- 
erate. This leads to an FBP for ()A.2p with the unknown domain of singularity, {(x, y) : n(x, y, t) =0}, 
and speciflc regularity properties; see [18^ I19| and references therein. Alternatively, finite-time forma- 
tion of non-smooth free boundaries with flat parts is a typical phenomenon for blow-up solutions of 
the reaction-diffusion PDEs (see the first equation in (jl.6p below) via extended semigroup theory. In 
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ID, optimal regularity of such C^'^-interfaces is well understood [34:\ Ch. 5]. For > 1, the regular- 
ity problem remains essentially open; see some estimates and examples in [34, p. 151]. As a formal 
extension, note that a nontrivial {u{x, t) ^ 0) proper convex solution exists for such PDEs with an 
arbitrarily strong (as u — > 0) absorption term, e.g., 

_ detp-^u _ lju 
"* ~ (l + |V«|2)3/2 « ' 

where a similar FBP occurs. Therefore, the parabolic operator of the Gauss curvature flow is extremely 
powerful, in the sense that it prevents a complete extinction (i.e., u{x,t) = for arbitrarily small 
t > 0; this can happen for many other parabolic PDEs). For any initial data with flat sides, {{x,y) : 
uo{x,y) = 0} 7^ 0, this proper solution of the FBP can be constructed by regular approximations of 
the equations and initial data by replacing e^/" i-^ min{^, e^^"}, with e > 0, (a uniformly Lipschitz 
continuous approximation), and uq ^ uq + e. Uniform a priori estimates for {u^} are obtained by 
local (near the interface) comparison with ID TW solutions or other radial sub- and super-solutions, 
[Ml Ch. 7]. 

In connection with other PDEs, let us also mention an unusual hyperbolic M-A equation 

UUtt + UttUxx - {Uxtf + 1 = 0, 

which occurs in ID hyperbolic mean curvature flows for closed curves Ftt = k'N — Vp on x [0, T), 
Vp = {Fts,Ft)T, where N and T are the unit normal and tangent vectors [56] (then u stands for the 
support function of F, the flow supports convexity, and finite-time shrinking to a point occurs). 
Concerning other nonlinearities, the elliptic M-A equation 

(A. 3) (detD'^u)^ = u <0, D^u> in n, -u = on (90, 

where is a bounded convex domain in M^, was derived by Loewner-Nirenberg [66] in the study of 
the metric of the form — ^ D^u {u is then treated as a section of a certain line bundle), and was proved 
to admit a unique solution; see [11] and earlier references therein. The general Hessian equation 
has the form 

(A.4) SkiD'^u) = {-uf, u<0 in n, u = on dfl, 

where is a baU in M^, A^ > 3, and Sk is given by the elementary symmetric function 

with {Aj} being the eigenvalues of the Hessian D'^u (so k = 1 and k = n correspond to the Laplace 
and the M-A operators, respectively). (|A.4p is known to exhibit the critical exponents 

— {N-2k)+ ' 

such that no smooth solution u < exists for p > 7(fe), and a negative radial solution exists for 
p E (0,7(/c)); see ^STj (nonexistence is proved by a Pohozaev-type inequality) and [12] for extensions. 
For k = 1, 7(1) = is the critical Sobolev exponent. The nonexistence result for the elliptic 

equation 

Au + uP = 

is associated with Pohozaev's classic inequality [76]. The exponents ^{k) above are to be compared 
with the critical ones 

-f{k) = j^^^ for elliptic PDEs - {-/\fu + \u\'p~^u = 0, 
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where the existence- nonexistence results are proved by higher-order Pohozaev's inequahties [77] apphed 
to general quasilinear 2A;th-order PDEs. For the M-A equation in a convex $7 C M.^ , 

det D'^u = p{x)g{—u) in fi, = on dQ, 

with positive nonincreasing g and p £ C°° also positive, a criterion of existence of convex negative 
solution is established in [TOl [7l] . Removable singularity theory for elliptic equations (these include 
M-A, Hessian, and Weingarten ones) 

F{D'^u) = i;{x,u,Du), 

where e.g., F{r) > F{s) + F{r — s) for r — s > (say, F is continuous and concave in a set oi N x N 
real symmetric matrices) and F{r) > 7(detr)^/^, is well developed; see [83] for further references and 
results. Note that a first removability theorem for the equation ()A.9|1 in is due to Jorgens (1955) 
|52j . See the beginning of Section |4] for an extra discussion of regularity issues for the M-A equations. 

As another standard nowadays direction of elliptic theory, let us mention recent research on existence 
of solutions of real (complex) M-A equations with infinite values (blow-up) on the boundary of a 
strictly convex bounded smooth domain 17 C M^: 

(A. 5) det D'^u{x) = g{x)f{u{x)) in Q, u{x) —>■ oo as x —>■ d^l; 

see [SOKMIES] for history, references, and results. Note that, for the semilinear elliptic equations 

(A.6) An = /(n), 

first results were obtained for f{u) = e" by Bieberbach in 1916 in 2D [4j and by Rademacher in 1943 
in 3D followed by the research for general nonlinearities f{u) by J.B. Keller and Osserman in 1957; 
see more historical details in [69], and Veron [89^ I90j and Labutin [60j for a more complete overview 
of other results. 

In particular and for instance, the above presentation suggests generalized second-order M-A par- 
abolic flows (these equations are formulated for u{x,t) being convex or "almost convex"), 

(A.7) ut = {detD^r ± i-uY, 

with some exponents m > and p G M, that generate many open problems concerning local existence 
of convex solutions, free-boundary (degeneracy set) propagation, extinction, and blow-up singularity 
patterns, etc. Other interesting models occur by choosing the elliptic operator as in (jX2]) . In radial 
setting, where (IA.7|) reduces to a ID quasilinear parabolic PDE, the interface equations and their 
regularity, moduli of continuity of proper solutions, waiting time phenomena, etc., are characterized 
by Sturmian intersection comparison techniques, [34^ Ch. 7]. For > 1, the majority of the problems 
are open, and particular exact solutions might be key. As we have shown, such models can be 
considered as natural counterparts of the PME with reaction/absorption, and of thin film (or Cahn- 
Hilliard-type, n = 0) models, such as (11. Oh . On basic properties of hyperbolic M-A equations, see [84^ 
Ch. 16] . Various parabolic and hyperbolic M-A equations admitting exact solutions on linear invariant 
subspaces are described in [37, Ch. 6]. 

Concerning blow-up or extinction behaviour in M-A flows that are main subject of the present paper 
and which often are not well-understood, as a simple illustration, a few types of singularity formation 
phenomena occurring on linear subspaces admitted by such Monge-Ampere operators are shown in 
[371 § 6.5]. As for (jl.5p . these are associated with invariant subspaces of the principal M-A operator 

(A.8) F2[u] = detD'^u = UxxUyy — (uxy)'^ in IK^- 
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Notice that quadratic polynomials p{x) occur in the celebrated result of the theory of elliptic M-A 
PDEs, establishing that any convex solution of the elliptic M-A equation 

(A.9) detD^u = 1 in 

is u{x) = p{x). This result is due to Jorgens (1954) for N = 2, Calabi (1958) for iV = 3, 4, and 5, 
and to Pogorelov (1978) for any N > 2 (see also |11) for a more general result). The same conclusion 
holds for the Hessian quotient equation 

Sk{D'^u) = 1 for u{x) < A{1 + |xp) strictly convex, 

with any 1 < < [2\. Similarly, if u{x,t) is a smooth solution of the parabolic PDE 

-ut detD'^u = 1 in X R+ , 

where u is convex in x, nonincreasing with t, and ut is bounded away from and — cxd, then u{x,t) = 
Ct + p(x); see [32] for the results and a survey. 

On higher-order M-A flows and blow-up. Fourth and higher-order M-A PDEs have been less well 
studied, though some of the equations correspond to classical geometric problems, and several general 
results have been established. We refer to [91j, where existence was established for the following class 
of mth-order fully nonlinear equations of the M-A-type in X Ml (m > 2) 

(A.IO) Ut = afj_^i^D^uD^u + (lower-order terms). 

Here fi and 1/ are multi-indices, and the matrix ||ci^|;/|| satisfies a positivity-type assumption for local 
existence. In |91] . a Riemann-Hilbert factorization condition appeared. In [85], '^-regularity esti- 
mates for fourth-order M-A equations were derived and an analogy of the Jorgens-Calabi-Pogorelov 
result for such PDEs was established. In [62], homogeneous fourth-order PDEs for affine maximal 
hypersurfaces were studied. Further references can be found in these papers. 

In order to formulate fourth-order Hessian equations in M?, let us write down the fourth differential 
of a C^-function u = u{x, y) quartic form 

d'^u = Uxxxxdx'^ + 4uj:j:j:ydx^dy + Guxxyydx'^dy'^ + 4uj;yyydxdy^ + Uyyyydy'^. 

This gives the catalecticant determinant (15. 2p . 

Using the operator F^lu] = detW^u, it is easy to construct some formal exact geometric flows. 
Clearly, F4 preserves the subspace of fourth-degree polynomials W7 = C{1, x'^,xy, y^, x^, x^y^, y^} and 
i<4 : Wj -^{1}- Therefore, the flow ut = F4[u] is global on Wj. The basic invariant subspace of 
sixth-degree polynomials is W = C{x°'y^ , < a + /? < 6}, on which blow-up may happen via a 
cubic DS. Singular patterns also exist for other fourth-order M-A-type models that are constructed in 
accordance to their second-order counterparts. Let us present two other examples with different types 
of evolution singularities. 

Let us discuss simple examples on extinction and blow-up for higher-order M-A flows. For instance, 
the M-A equation (cf. (jA.ip with the known extinction behaviour) 

^t = -ditS^ inM2xM+ 

admits solutions on the subspace W4, u{x,y,t) = Ci{t) + C2{t)x'^ + Cs{t)x'^y'^ + C4{t)y^ . Then \D^u\ = 
64(36C2C3C4 — C3). The corresponding fourth-order DS yields solutions in separate variables 

(A.ll) u{x,y,t) = Ci{t){l + Ax^ + Bx^y^ + Cy'^) , 
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where A, B, and C are positive constants satisfying 7 = 64:{3GABC — B"^) > by the convexity 
assumption on initial data. Here Ci solves the ODE C[ = This gives finite-time extinction 

with the rate 

= [l(r-t)]^ ^0 as t^T~. 

Vice versa, 

ut = uVdetD^n in M? x M+ 
admits solutions (jA.lip driven by the ODE C| = C^^"^ with blow-up, 



Similar singularity phenomena are traced for the corresponding hyperbolic M-A flows on these sub- 
spaces. A number of typical conclusions for the second-order Hessian flows can be extended to this 
fourth-order, as well as higher-order, though the well-posedness of such parabolic or hyperbolic PDEs 
in classes of "convex" functions is a difficult open problem. 

Finally, notice that the homogeneous equation = is a direct sum of two identical copies of 

the second-order M-A equation (see [29j ) 



(A.12) 



VxxVyy 



0. 



Similarly, the sixth-order equation \D^u\ = with the operator 



(A.13) 



detD^u = det 



u. 



x'y- 



■U60 ^^51 '"42 ^^33 
■"51 ^^42 -"33 U2i 
U42 U33 U24 Ul5 
U33 U24 Ui^ U06_ 

decouples into three copies of ()A.12p . |29| . Possibly, this means that some problems with such higher- 
order M-A operators are associated with the second-order ones. In particular, the inhomogeneous 
equations \D'^u\ = 1, \D^u\ = 1 might be handled by reduction to second-order equations, and a result 
associated with the Jorgens-Calabi-Pogorelov theorem might be expected (though some basics of such 
PDEs remain obscure). 

Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, UK 
E-mail address: cjb@maths.bath.ac.uk 

Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, UK 
E-mail address : vagSmaths . bath . ac . uk 



43 



